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Identified 7T , K ' , p and p transverse momentum spectra at mid-rapidity in 
^/sjvjv = 130 GeV Au-Au collisions were measured by the PHENIX experiment at RHIC as a 
function of collision centrality. Average transverse momenta increase with the number of participat- 
ing nucleons in a similar way for all particle species. Within errors, all mid-rapidity particle yields 
per participant are found to be increasing with the number of participating nucleons. There is an 
indication that K + ^~ , p and p yields per participant increase faster than the 7r +//_ yields. In central 
collisions at high transverse momenta (pt ~ 2 GeV/c), p and p yields are comparable to the 
yields. 

PACS numbers: 25.75.Dw 



We report first results on identified 7r+/~ , K + /~, p 
and p production as a function of collision centrality in 
y?SNN — 130 GeV Au+Au collisions, measured by the 
PHENIX experiment at the Relativistic Heavy Ion Col- 
lider (RHIC). The PHENIX objective is to search for 
signatures of deconfincment and chiral symmetry restora- 
tion and to study the transition from normal to decon- 
fined nuclear matter by utilizing a wide variety of probes. 

Early RHIC results show that the transverse energy 
density and particle multiplicities are considerably higher 
than previously observed in relativistic heavy-ion colli- 



sions |l||2|. Measured energy densities extend into the 
region predicted to be favorable for the formation of a 
quark-gluon plasma Identified hadron spectroscopy 
provides a tool for studying reaction dynamics beyond 
that of global event characterization. The yields of 
hadrons reflect the particle production mechanism, while 
spectral shapes are sensitive to the dynamical evolution 
of the system. The mass and centrality dependence of 
the spectra can help differentiate between competing the- 
oretical descriptions such as collective hydrodynamical 
expansion or transverse momentum (pt) broaden- 
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ing in the partonic stage of the reaction Q. Addition- 
ally the relative yields of baryons and mesons at high 
(j>t — 2 GeV/c) may give insight into baryon number 
transport j?j and the interplay between soft and hard 
processes. 

The PHENIX detector has diverse particle identifica- 
tion (PID) capabilities 0], including excellent hadron 
identification over a broad momentum range. This mea- 
surement was performed using a portion of the east 
central-arm spectrometer, covering pseudo-rapidity \rj\ < 
0.35 and A(f> = 7r/4 in azimuthal angle. 

The collision z-vertex and the timing system's start 
signal are generated by the Beam-Beam Counters (BBC); 
two arrays of quartz Cherenkov radiators which surround 
the beam axis covering r\ = ±(3.0 — 3.9). The track- 
ing system includes a multi-layer focusing drift cham- 
ber located outside an axially-symmetric magnetic field 
at a radial distance between 2.0 m and 2.4 m followed 
by a multi-wire proportional chamber with pixel-pad 
readout (PCI) |J. Pattern recognition in the DC is 
based on a combinatorial Hough transform in the track 
bend plane [^0| . The polar angle of the track is deter- 
mined by PCI and the collision z-vertex. A track model 
based on a field-integral look-up table determines the 
charged particle momentum and the path length to the 
time-of-flight (TOF) wall. The momentum resolution is 
Sp/p ~ 0.6% ffi 3.6% p (GeV/c). 

The timing system stop signal for each particle is mea- 
sured by the TOF scintillator wall, located at a radial 
distance of 5.06 m, resulting in a flight-time measure- 
ment with a resolution of a ~ 115 ps. Reconstructed 
tracks are projected to the TOF and matched with hits in 
the scintillator slats using a momentum-dependent search 
window determined by multiple scattering and the mo- 
mentum resolution. A velocity dependent energy loss 
cut based on a Bethe-Bloch parameterization is applied 
to the measured TOF pulse height. Combining the mo- 
mentum and flight-time, we reconstruct the particle mass 
and select particles by applying 2a momentum dependent 
cuts in mass-squared. 

Corrections for geometrical acceptance, decay- in- flight, 
momentum resolution and reconstruction efficiency are 
determined using a single-particle full GEANT Monte 
Carlo (MC) simulation. The acceptance correction as- 
sumes that the spectra are flat in azimuth and in rapid- 
ity for \y\ < 0.5 [yj. Fiducial area cuts, energy loss and 
hit-track matching cuts are applied consistently in sim- 
ulation and data. In peripheral events the track recon- 
struction efficiency is ps 98%. As the centrality increases, 
the efficiency is reduced due to the increased detector oc- 
cupancy. Multiplicity dependent corrections are obtained 
by embedding simulated tracks into real events. Track- 
by-track corrections are applied taking into account the 
event centrality and the particle species. The efficiency 
for 7r + / - in the most central events is ps (68 ± 6)%, 
independent of momentum. In the case of overlapping 



hits in the TOF wall, the earliest pulse reaching each 
photo-multiplier is recorded. This favors the faster par- 
ticles; hence, in central events heavier particles suffer an 
additional reconstruction inefficiency of ps 4%. Correc- 
tions for feed-down from weak decays are not applied. A 
MC simulation is used to estimate the probability for 
reconstructing protons from A decays as prompt pro- 
tons. Within the PHENIX acceptance this probability 
is w 50% at p T = 0.5GeV/c, ps 32% at p T = 1 GeV/c 
and ps 12% at px > 2 GeV/c. Taking ^ = 1 as an upper 
limit, we estimate 33% as the upper limit of weak decay 
contribution to the reported p, p yields and maximal 16% 
change of measured <px>- 

About 140,000 minimum bias events, representing 92 ± 
4% of the total inelastic cross-section of 6.8 b were 
analyzed. This sample was subdivided into five central- 
ity classes: - 5%, 5 - 15%, 15 - 30%, 30 - 60% and 
60 — 92%, using the BBC and Zero-Degree-Calorimeters 
for event characterization jjj . For each class, the average 
number of nucleons participating in the collision (N par t) 
is obtained from a Glauber model calculation [Q . 

Fig. [l] shows the invariant yield as a function of pt 
for 7r + , K + ,p (left panel) and ir~,K~,p for three cen- 
trality selections. Error bars in the figure are the com- 
bined statistical errors in the data and the corrections. 
Systematic uncertainties from acceptance, multiplicity- 
dependent efficiency corrections and PID cuts result in 
an overall systematic error in the absolute normalization 
of ps 11% for all species. As a consistency check, we 
have added all identified charged hadron spectra in the 
Pt region where PID is available for all species and com- 
pared to the PHENIX charged hadron measurement jl2) . 
The results agree to better than 10% for all centralities. 
Additionally, the PHENIX tt° spectra flf] and the 
spectra in the region of overlapping px are within ps 10% 
for the central and ps 25% for the peripheral selections, 
which is within the systematic errors of the two measure- 
ments, p results presented here are also in agreement 
(within errors) with recent publication |13| ]. 

In peripheral events the 7r + /~ spectra exhibit a concave 
shape, well described by a power-law parameterization as 
observed in hadron-hadron collisions E3. With increas- 
ing centrality the curvature in the spectra decreases, lead- 
ing to an almost exponential dependence on px for the 
most central events. Over the measured px range, the 
K + /~ spectra can be described by an exponential dis- 
tribution either in px or in rax — yPx + m2 j while the 
p and p spectra can be described either as a Boltzmann 
or an exponential distribution in mx- The slopes of the 
mx spectra flatten and the mean transverse momentum 
(<px>) increases with particle mass and with central- 
ity. This behavior has been previously observed in lower 



energy heavy ion collisions at the BNL-AGS |15| and at 
the CERN-SPS ^§0] and was attributed to collective 
radial motion (flow). 
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At lower energies, it is not uncommon for the proton 
yields to equal or exceed the ir + yields, since many of 
the protons come from the initial state. A new feature 
observed for the first time at RHIC is that in central 
collisions at px ~ 2GeV/c the p yields are compara- 
ble to the 7T~ yields. Positive and negative hadrons 
behave in a similar way. In central events the proton 
yields approach the n + spectra around px — 1.6GeV/c. 
As the centrality decreases, this happens at larger px- 
In peripheral events, p and p spectra are below the 
7r + /~ spectra over the whole measured px range. Since 
anti-protons are not as numerous as the protons, p and 
7r~ yields become comparable only at the high end of 
the measured pion px range in the most central colli- 
sions. We note that in pp Q and pp @ collisions at 
y/s = 23-63GeV and yfi = 300- 1800 GeV respectively, 
the p/tt~ ratio steadily increases with px up to ~ 0.33 
at pr — 1.5GeV/c nearly independent of y/s. Data on 
baryon/meson ratios at higher px are only available from 
pp collisions at ISR energies (y/s = 23 — 63 GeV) and 
show that above px — 1.5GeV/c the p/^~ ratio rises to 
~ 0.4 at px — 2GeV/c and then drops, as expected if 
valence quark jet fragmentation is the dominant produc- 
tion mechanism for high px hadrons. In central collisions 
at RHIC the high px (anti)proton/pion ratios are ~ 1- 
much larger than in pp collisions. 

Hydrodynamic calcualtions with fixed freeze-out tem- 
perature B or with freeze-out modeled using a hadronic 
cascade [[| suggest that hydrodynamic expansion is re- 
sponsible for the baryon dominance at high px- How- 
ever, protons/anti-protons produced via a baryon junc- 
tion mechanism combined with jet-quenching in the pion 
channel are shown to exhibit the same effect jjj. In- 
trinsic px broadening in the partonic phase caused by 
gluon saturation j^] gives yet another alternative expla- 
nation. The above models have similar predictions in the 
px range measured here, but show different behavior at 
higher px- New data with broader px range is needed in 
order to distinguish between currently available theories. 

To quantify the centrality and mass dependencies of 
hadron production, we determine the mean (<px> shown 
in Fig. ||) and the integral (dN/dy shown in Fig. ||) of 
the px distributions. Both quantities require extrapola- 
tion of the spectra below and above the measured range. 
For each particle species, we use at least two functional 
forms consistent with the data, as outlined above and the 
results presented are averaged between two fits. The frac- 
tion of the yield in the extrapolated regions is estimated 
30 ± 6% for 7T+/-, 40 ± 8% for K + /~ and 25 ± 7.5% for 
p/p- Combining systematic uncertainties in the extrap- 
olation fractions and the estimated background under 
mass-squared peaks (2%, 5% and 3%) yield systematic 
uncertainties in the measured <px> of 7%, 10% and 8% 
for 7r + / - , K + /~ and p/p, respectively. Further combin- 
ing the above uncertainties, which only affect the shape 
of the spectra, with the 11% uncertainty in the abso- 



lute normalization gives systematic uncertainties of 13%, 
15% and 14% in the measured dN/dy for 7r +/ ~, K +/ ~ 
and p/p, respectively. We note that after converting to 
dN/d-q, the sum of the identified charged hadron yields 
agrees with the previously published PHENIX results on 
total charged multiplicity dN c h/dr] 111 within 5%. 

Fig. U shows the <px> as a function of N part for 
-K + ,K + ,p (left panel) and Tr~,K~,p. Filled points are 
this measurement; open points at N par t — 2 are interpo- 
lations to y/s = 130 GeV obtained from pp and pp data 
at lower and higher energies [^9|, respectively. In 
peripheral Au-Au collisions at RHIC 7r + / - and K + /~ 
exhibit similar <px> to those in pp collisions, but pro- 
tons and anti-protons have significantly higher <px>, in- 
dicating that nuclear effects are important even at small 
Npart- For all measured particles species, <px> increases 
by as 12 — 14% from the first to the second centrality bin 
(i.e. N part 14 to 79). Above Npart = 100 the ir +/ ~ and 
K + /~ <px> appear to saturate, whereas the p and p 
<px> rises slowly. However, we note that going from 
peripheral to the most central event class the overall in- 
crease in < px > is ~ 20 ± 5% independent of particle 
species. 

Fig. |^ shows the yields per participant versus N par t- 
Error bars include statistical and multiplicity dependent 
systematic errors. Systematic uncertainties in N par t that 
can move all curves independent of mass are shown with 
bands around the positive hadron points. Total system- 
atic uncertainties in N par t and the yields at each N part 
are listed in Table [jj. 

For all particle species, the yield per participant in- 
creases with Np ar t. As for <px>, most of the increase 
occurs between the two most peripheral centrality selec- 
tions (i.e. Npart 14 to 79). However, in contrast with 
the centrality dependence of <px>, we see an indication 
that the total increase in the yields per participant dif- 
fers among particle species as we go from peripheral to 
central events. The pion yield per participant rises by 
21% ± 6%(stat) ± 8%(syst). The kaon yields per par- 
ticipant rise faster: 94% ± ll%(siai) ± 2Q%(syst) and 
66% ± 12%(stat) ± 20%(syst), for K+ and K~ , respec- 
tively. Similar trends in the centrality dependence of 
strangeness production have been observed in lower en- 
ergy heavy ion collisions at the BNL-AGS J2(]] and at the 
CERN-SPS 0. It is interesting to note that at RHIC p 
and p yields per participant behave similarly to the K + /~ 
yields and also rise faster than the pions with increasing 
Npart- The increase is 58% ± 5%(stat) ± 16%(syst) and 
72% ± 9%(stat) ± 20%(syst), respectively. The similar 
centrality dependence in p and p yields per participant 
indicates that baryon/anti-baryon pair production is the 
dominant source of protons and anti-protons alike. 

In heavy ion collisions at AGS energies |^| the p pro- 
duction is close to threshold, the yields per participant 
are lower than in pp collisions and decrease from periph- 
eral to central collisions, probably due to annihilation. 
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At the SPS, the p yield per participant is larger than the 
pp value and has almost no centrality dependence p3| . 
At RHIC, the total yield of anti-protons at mid-rapidity 
in central Au-Au collisions is a factor of ~ 1000 larger 
than at the AGS (2^] and nearly an order of magnitude 
above that in Pb+Pb collisions at CERN Q. Most of 
the increase is due to the y/s dependence of baryon/anti- 
baryon pair production, however the yield per participant 
rises noticably from peripheral to central collisions. 

In conclusion, an intriguing new behavior in identified 
hadron production at RHIC is reported. In central Au- 
Au collisions the anti-protons yield is comparable to the 
7r~ at high pt - a behavior never observed before in el- 
ementary or in heavy-ion collisions. <pr> rises with 
centrality similarly for all particle species, while K + /~, 
p and p yields per participant increase somewhat faster 
than the 7r + /~ yields. 
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FIG. 1. Transverse momentum spectra measured at 
mid-rapidity for n + ,K + ,p (left) and tt~ , K~ ,p at the three 
different centrality selections indicated in each panel. The 
symbols indicated in the top panels apply for all centrality 
selections. 
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FIG. 2. Average transverse momentum for n + , K + ,p (left) 
and 7r~, K~ ,p as a function of the number of nucleons par- 
ticipating in the collision N par t- The error bars represent the 
statistical errors. The systematic errors are discussed in the 
text. The open points are interpolations from pp and pp data, 
see text for details. 
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FIG. 3. dN I 'dy\ y=0 per participant for ir + ,K + ,p (left) and 
ir~ , K~ ,p as a function of N par t- The error bars include sta- 
tistical and systematic errors in dN/dy. The dashed lines 
around the positive hadrons show the effect of the systematic 
error on N part which affects all curves in the same way. 



TABLE I. Integrated hadron ( ty ± ,K ± , p and p ) yields 
at mid-rapidity for five centrality classes (see text) identi- 
fied by the indicated number of participants N par t- The er- 
rors on N par t are the systematic errors. The errors listed for 
<dN / dy\ y = > are statistical. The systematic errors are 13%, 
15% and 14% for n ± , K ± , p and p, respectively. 

N part 348+10.0 271+8.4 180+6.6 79+4.6 14+3.3 

276+3 216+2 141+1.5 57+0.6 9.6+0.2 

7T" 270+3 200+2 129+1.4 53.3+0.6 8.6+0.2 

K + 46.7+1.5 35.0+1.3 22.2+0.8 8.3+0.3 0.97+0.11 

K~ 40.5+2.3 30.4+1.4 15.5+ 0.7 6.2+0.3 0.98+0.1 

p 28.7+0.9 21.6+0.6 13.2+0.4 5.0+0.2 0.73+0.06 

p 20.1+1.0 13.8+0.6 9.2+0.4 3.6+0.1 0.47+0.05 
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